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We study multifractal spectra of critical wave functions at the integer quantum Hall plateau 
transition using the Chalker-Coddington network model. Our numerical results provide important 
new constraints which any critical theory for the transition will have to satisfy. We find a 
non-parabolic multifractal spectrum and determine the ratio of boundary to bulk multifractal 
exponents. Our results rule out an exactly parabolic spectrum that has been the centerpiece in 
a number of proposals for critical field theories of the transition. In addition, we demonstrate 
analytically exact parabolicity of related boundary spectra in the two-dimensional chiral orthogonal 
'Gade-Wegner' symmetry class. 
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The physics of the quantum Hall effect has been an ex- 
citing area of research for more than two decades |j], |) . 
While much progress has been made in this area, the 
identification of an analytically tractable theory describ- 
ing the critical properties at the transitions between the 
plateaus in the integer quantum Hall (IQH) effect has 
been elusive ever since p|. These quantum phase tran- 
sitions are famous examples of (Anderson) localization- 
delocalization (LD) transitions driven by disorder. The 
diverging localization length plays the role of a correla- 
tion length in non-random continuous phase transitions, 
known to be described by conformal field theories in two 
dimensions (2D). It is natural to expect that effective 
(field) theories describing IQH plateau transitions should 
generally also possess conformal symmetry (cf. M). 

Many attempts have been made in the past to identify 
an analytically tractable description of the IQH plateau 
transition and, more recently, Wess-Zumino (WZ) field 
theories defined on a certain supermanifold were conjec- 
tured to provide such a description p, |fj, In. (Similar 
theories have also appeared in the context of string prop- 
agation in Anti-de Sitter space-time ||.) These propos- 
als focussed solely on bulk obscrvables, i.e., on physical 
quantities measured in a sample without any boundaries. 
In this Letter, we provide important new constraints that 
arise when one studies the scaling behavior of wave func- 
tions near the boundaries of a sample. Any proposed 
candidate theory for the plateau transitions will have to 
be consistent with our numerical results for the boundary 
multifractal spectrum. 

At LD transitions, critical wave functions obey scale- 
invariant, multifractal (MF) statistics, namely, disorder- 
averaged moments of wave functions have a power-law 
dependence on the linear dimension L of the system || : 

|2\9 
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derlying field theory flC)|j , can be defined for points r in 
the bulk (x = b) of the sample, A* or near its boundary 
('surface': x = s) @, ||], A". The prefactor C x q {V) in 
Eq. (|l|) depends on q and, in general, on L if we include 
the possibility of corrections to scaling. Both sets of MF 
exponents satisfy the symmetry relation ||12| 



A x = Af. 



(2) 



(in some interval Q around q = 1/2). 

Equivalently, the MF wave functions can be character- 
ized by the so-called singularity spectra f x (a x ) related 
to A x by a Legendre transform: f x (a x ) = (a x — 2)q — 
A x + d x , a x - 2 = dA x /dq, and d b = 2, d s = 1. The 
expon ent a x de scribes the scaling of typical wave func- 
tions: In ^(r)! 2 ~ — a x InL, as can be seen by taking the 
q derivative in Eq. (fil) at q = 0. 

Work emerging jljj [LJ] from Ref. || led to the con- 
jecture that the proposed theory would give rise to an 
exactly parabolic bulk MF spectrum for the IQH transi- 
tion 



l"q(l-q), 



(3) 



The MF exponents A x which are related to ('anoma- 
lous') scaling dimensions of certain operators in an un- 



reminisccnt of analytically obtained MF spectra for 
Dirac fcrmions in, e.g., random abclian gauge potentials 
|15| , |lq] . In those models the parabolicity of the MF spec- 
trum can be understood through a reformulation of the 
problem in terms of free fields. 

Previous numerical studies |l3[ of wave function statis- 
tics at the IQH transition appeared to exhibit a bulk 
MF spectrum that was indeed well described (with an 
accuracy of ~1%) by a parabolic fit (||) with 7 6 = 
0.262±0.003, seemingly providing support for the conjec- 
tures advanced in Ref. [|| ||, ItJ . (In Ref. Q the results 
are presented in terms of f (a). For a parabolic MF 
spectrum (0), f b (a b ) is also parabolic, with a maximum 
at a b = 7 b + 2.) 



Besides its conjectured relevance || to the IQH tran- 
sition, the above-mentioned WZ theory is known to de- 
scribe transport properties of a disordered electronic sys- 
tem in a different universality class Jl7|, |lq] (the chiral 



unitary 'Gade-Wegner' class AIII of 20, til) which pos- 
sesses an additional discrete (chiral) symmetry [E0( , not 
present in microscopic models for the IQH transition. 
Well-known microscopic realizations of field theories in 
class AIII are random bipartite hopping models, and cer- 
tain network models |0, O, 19[. The theory possesses 
a line of fixed points, with continuously varying critical 
properties parametrized by the critical longitudinal DC 
conductivity. (It was argued in Rcf. || that for a partic- 
ular value of this continuous parameter the WZ theory 
would provide a description of the IQH transition.) 

In this paper we obtain two kinds of results. First, we 
provide results of extensive numerical work on the MF 
exponents at the IQH transition both at a boundary (A;?) 
and in the bulk (A*). Based on these numerical results 
quadratic behavior in q is ruled out for both quantities. 
Deviations from the parabolic form (pf) are found to be 
much larger in the MF exponents A* at a boundary. Here 
it is important to note that in complete analogy to the 
bulk, the above conjectures would also yield a quadratic 
dependence on q of the boundary MF exponents A^. We 
further determine the ratio A*/A^ over a range of q. 
Accounting for this ratio is an important constraint on 
any proposed critical theory for the transition. 

Secondly, we demonstrate analytically the exact 
parabolicity of boundary spectra, not for the chiral uni- 
tary class AIII, but for the related time-reversal invariant 
version, the chiral orthogonal 'Gade-Wegner' class BDI 
jl7], [l8| [ll], ||(J . We expect such parabolicity to also hold 
in the chiral unitary symmetry class. 

We begin with the numerical part. Here, we study the 
multifractality of critical wave functions in a way similar 
to Ref. jl3| , with the goal of numerically determining the 
rescaled anomalous exponents 



7j = A|Ml-g), 



(4) 



both for x = s (boundary) and x = b (bulk). 

For the case of boundary exponents we consider the 
critical Chalker-Coddington network model (CCNM) |2§] 
with 4L 2 links placed on a cylinder. The dynamics of 
wave functions on links of the network is governed by 
a unitary evolution operator U. For each disorder re- 
alization, we numerically diagonalizc U and retain one 
critical wave function whose eigenvalue is closest to 1. 
The largest system size we studied was L = 180, and 
the ensemble average was taken over 3 x 10 5 samples for 
L = 50,60, 5 x 10 5 samples for L = 80, and 2 x 10 5 
samples for L = 120, 180. 

We obtain the anomalous dimensions A* from Eq. 
(H). The boundary wave function coarse-grained over 
each plaquette along the boundary, is substituted into 
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FIG. 1: (Color online) System size dependence of 

Dg(L)/q(l — q) defined in the main text for q — —0.2 (■), 
0.2 (•), 0.5 (a), 0.8 (o), and 1.2 (□). 7* is calculated by 
linear fitting taking into account only for larger system sizes 
(L ^ 50) indicated by the vertical dashed line between L = 40 
and 50. We also show A s — 2 defined in the main text (x). 

the left-hand side of Eq. (||), where the over line de- 
notes ensemble and spatial averages along the boundary. 
Taking the logarithm, we numerically obtain D S (L) = 

{q\n\i>{r)\ 2 - In \ip{r)\ 2 i) / lni = A* - lnC*(L)/lnZ, 
and plot this quantity as a function of 1/lnL in Fig. 
[y. We see that corrections to scaling are significant for 
small systems (L < 30). Therefore, we used our numeri- 
cal data for L ^ 50 only to extract 7;? by linear fitting. 

Independently, wc numerically obtain a q and f s from 

\^{r)\ 2( i\-a\Hr)\ 2 /\ij:(r)\ 2 <i a s q In L, and ln^r)! 2 ? ~ 

\f s (a q ) — a q q — ds\ln.L, using our numerical data for 
L ^ 50. For example, the expone nt ag is obtained 
by linear fitting (Fig. |lj), A s = -In \ip{r)\ 2 /lnL = 
ag + const./ lnX, which yields a s = 2.386 ± 0.004. 

We show in Fig. g(a) the rescaled boundary anoma- 
lous dimension 7;? (red filled circles) obtained from this 
analysis. We see clearly that 7! is not constant, imply- 
ing that the boundary MF spectrum A* is not parabolic. 
The change in 7* over the interval < q ^ 1/2 is about 
4 ~ 5% and is significantly larger than the error bars. 
This provides the strongest numerical evidence against 
the parabolicity of the MF exponents. 

Shown in the same figure by blue open circles is the 
mirror image of 7* with respect to q = 1/2, 7f_„. We see 
that the symmetry relation (0) is satisfied within error 
bars for < q < 1. The rescaled anomalous dimension 
7* approaches ag — 2 (the horizontal line) at q = 0, 1, 
indicating that the two independent calculations of ag 
and A* are consistent. 

We have also computed the bulk anomalous dimension 
A* using the CCNM on a torus. In this case the over- 
line in Eq. ([j]) implies both the ensemble and the spatial 
average over the whole torus. Wave functions are coarse- 
grained on each plaquette. We have employed the same 
fitting procedure as in the boundary case. The biggest 
system size we examined for the bulk analysis is L = 270. 
The number of samples over which we took the average 
is 5 x 10 5 for L = 50, 3 x 10 5 for L = 60, 80, 2 x 10 5 for 
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FIG. 2: (Color online) (a) Rescaled boundary MF exponents 
7| (•) and 7l _ 9 (o). The curve is 0.370 + 0.042(g - 1/2) 2 , 
obtained by fitting the data for 7° in < q < 1 to a parabolic 
form. The horizontal solid line shows oiq — 2 = 0.386 ± 0.004 
with error bars indicated by dashed lines, which is consistent 
with liniq^o,! 7,. (b) Rescaled bulk MF exponents 7J (•) and 
7i_ q (o). The curve is 0.2599 + 0.0065(g - 1/2) 2 obtained 
by fitting the data for 7^ in < q < 1 to a parabolic form. 
The horizontal solid line shows a b - 2 = 0.2617 ± 0.0006 with 
error bars indicated by dashed lines, (c) Ratios Jq/jq (•) 
and 7i_,j/7i_q (o). As above, the curve is obtained from the 
parabolic fits for 7|' i> , which amounts to quartic approxima- 
tions for A s ' b . 



L = 120, 4 x 10 4 for L = 180, and 2 x 10 4 for L = 270. 

Figure 0(b) shows the exponents 7* together with 
their mirror image. The symmetry relation (0) is again 
satisfied for < q < 1 within error bars, which provides 
confirmation that our results are reliable. We see clearly 
that 7^ has q dependence, although it is weaker than that 
of 7^; compare the vertical scales of Fig. g(a) and (b). 

The ratio "fq/jq is shown in Fig. 0(c) and is seen to 
be clearly dependent on q. Any candidate theory for the 
IQH transition needs to be consistent with this ratio, and 
in particular its dependence on q. (Note that for a free 
field this ratio would be equal to 2, and independent of 

9 0, Ho 

Figure ||(a) shows a? as a function of q. The data sig- 
nificantly deviate from linear behavior that would result 
if A x were strictly parabolic (constant 7*). In Fig. |](b) 
we show the corresponding singularity spectra f x (a x ) as 
functions of q. (Combining the data from the two panels 
would result in f x (a x ) as functions of a x .) For q > 1.5 
where f s (a s ) < 0, the moments \ip(r)\ 2 i are dominated 
by rare events, and thus accurate numerical calculation of 




FIG. 3: (Color online) (a) a s q (•) and a q (o) as functions 
of q; (b) fq (•) and fq (o) as functions of q. The solid and 
dashed curves on both panels are obtained from the parabolic 
approximations to 7, (that is, quartic approximations to A^). 
Notice that q^ significantly deviate from straight lines which 
would follow from strictly parabolic A^ (or constant 7^). 

MF exponents becomes more difficult j9|. This explains 
the discrepancy between the (red) data points and the 
solid curves for q > 1.5 in Fig. |[ As /| > at q > — 1, 
we expect that our numerical results of f s {a s ) should be 
more reliable at q « —1 than at q ~ 1.5, as evidenced by 
the agreement between the red dots and the solid curve 
in Fig. g. The curve suggests termination of f s (a s ) H 
to occur at gw 2.2. 

We now proceed to present our second (analytical) 
result. We first recall that the theory conjectured in 
Rcf. Q to describe the IQH plateau transition is a WZ 
model with global psl(2|2) (super-) symmetry. This the- 
ory possesses two coupling constants, one denoted by 
/ for the kinetic term and another denoted by k for 
the WZ term (k = 1 in Ref. H), in standard notation 
P, E4fl. One can think of this theory as a perturba- 
tion of the rather well understood p4, p5| Kac-Moody 
(KM) point characterized by the condition f~ 2 = k, per- 
turbed by a term in the action of the form p4, Gm SS = 
(X/k 2 )Jd 2 z<j) ab (z,z)J a (z)J b (z) 1 where A = f~ 2 - k. 
Here <j) a b is the KM primary field in the adjoint repre- 
sentation of psl(2|2), J° and J b are the left/right chiral 
components of the psl(2|2) Noether currents [E4|, and A 
parametrizes the line of fixed points, mentioned above. 

The conjectured link Q between the WZ model and 
the IQH transition can be formulated through the notion 
of the point contact conductance (PCC) |27|. The PCC 
is a statistically fluctuating quantity, (i) On one hand, 
the scaling dimension X q of the g-th moment of the PCC 
at the IQH transition has been proven to be simply re- 
lated to the exponent A q Q Q, X q = 2A, (\q\ < 1/2). 
(ii) On the other hand, the scaling dimension x q of the 
operator in the WZ model carrying the same representa- 
tion of the global psl(2|2) symmetry as the g-th moment 
of the PCC in the CCNM (possessing the same psl(2|2) 
symmetry), was conjectured to be a quadratic func- 



tion of q. (iii) If one combines (i) and (ii), and if one 
assumes X q = x q (following the conjectured description 
of the IQH transition by the WZ model), then the wave 
function exponents A q at the IQH transition would be 
quadratic functions of q, as in Eq. (||). 

As already mentioned, this WZ theory is known to de- 
scribe transport properties of the chiral unitary class AIII 
J17], |TM |l9| |, lacking time-reversal symmetry. Below we 
demonstrate the correctness of the conjecture made in 
item (ii) above, at a boundary, and for the time-reversal 
invariant version of the AIII model, the chiral orthogonal 
class BDI 1 17, pC|. Just as its cousin with broken time- 
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reversal symmetry, the chiral orthogonal theory also pos- 
sesses a line of fixed points. Transport properties along 
this line can be described [|l8| by the perturbation of the 
KM point of the psl(2|2)-invariant WZ theory described 
above, when the field <p a b is replaced by the Kronecker 
delta, 4> a b — ► o~ab- Denote the corresponding coupling 
constant by At. Consider the theory in the upper half 
plane where the system simply ends at the boundary, and 
an operator of scaling dimension (a 'conformal weight') 
x s (Xt) on the boundary. At the KM point, where the per- 
turbation vanishes, At = 0, such an operator is described 
by a representation p of the global psl(2|2) symmetry. It 



(2) 



■,(2) 



is known |24j that x s p {\ t = 0) = Cf'/k, where C, 
the quadratic Casimir invariant in the representation p. 
It turns out to be straightforward (28[ to compute the 
change of the scaling dimension, order by order in the 
bulk coupling constant At, yielding a geometric series. 
The result is simply x s p (\ t ) = C { p 2) /(k + X t ). Note that 
for the (continuous series) representation p of psl(2|2) in 
which the g-th moment of the PCC at the IQH transi- 
tion transforms, one has C p = q{\ — q). This proves our 
claim that the spectrum of scaling dimensions x s — » x s 
of corresponding boundary operators in symmetry class 
BDI is a strictly quadratic function of q. 

In summary, our numerical results clearly demonstrate 
that both, the boundary and the bulk MF spectra, A* 
and A' significantly deviate from parabolicity, and that 
their g-dependent ratio is significantly different from 2. 
(These conclusions were recently also reached, indepen- 
dently, by Evers, Mildcnbcrger, and Mirlin p9[ .) These 
results for the bulk as well as the boundary MF spectra 
impose important constraints on any analytical theory 
for the IQH plateau transition. Furthermore, we have 
demonstrated analytically exact parabolicity of related 
boundary spectra in the 2D chiral orthogonal 'Gadc- 
Wegner' symmetry class BDI. 
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